arXiv: 1502.00914v2 [gr-qc] 9 Apr 2015 


Directed searches for continuous gravitational waves from binary systems: 
parameter-space metrics and optimal Scorpius X-1 sensitivity 

Paola Leaci^’^’El and Reinhard Prix^’0 

^Max-Planck-Institut fur Gravitationsphysik, Albert-Einstein-Institut, D-14476 Golm, Germany 
^Dip. di Fisica, Universitd di Roma “Sapienza”, P.le A. Moro, 2, 1-00185 Rome, Italy 
^ Max-Planck-Institut fiir Gravitationsphysik, Albert-Einstein-Institut, D-30167 Hannover, Germany 

(Dated: April 10, 2015) 

We derive simple analytic expressions for the (coherent and semi-coherent) phase metrics of 
continuous-wave sources in low-eccentricity binary systems, for the two regimes of long and short 
segments compared to the orbital period. The resulting expressions correct and extend previous re¬ 
sults found in the literature. We present results of extensive Monte-Carlo studies comparing metric 
mismatch predictions against the measured loss of detection statistic for binary parameter offsets. 

The agreement is generally found to be within ~ 10% — 30%. As an application of the metric 
template expressions, we estimate the optimal achievable sensitivity of an Einstein@Home directed 
search for Scorpius X-1, under the assumption of sufficiently small spin wandering. We find that 
such a search, using data from the upcoming advanced detectors, would be able to beat the torque- 
balance level mm up to a frequency of ~ 500 — 600 Hz, if orbital eccentricity is well-constrained, 
and up to a frequency of ~ 160 — 200 Hz for more conservative assumptions about the uncertainty 
on orbital eccentricity. 


I. INTRODUCTION 

Continuous gravitational waves (CWs) are a promis¬ 
ing class of signals for the second-generation detectors 
currently under construction: advanced LIGO (aLIGO) 
[3], advanced Virgo H] and KAGRAj^. These signals 
would be emitted by spinning neutron stars (NSs) subject 
to non-axisymmetric deformations, such as quadrupolar 
deformations (“mountains”), unstable oscillation modes 
(e.g. r-modes) or free precession. For a general review of 
GW sources and search methods, see for example [5]. 

A particularly interesting type of potential GW sources 
are NSs in low-mass X-ray binaries (LMXBs), with Scor¬ 
pius X-1 being its most prominent representative [7]. 
The accretion in these systems would be expected to 
have spun up the NSs to the maximal rotation rate of 
> 1000 Hz [5]. All the observations performed to date, 
however, show they only spin at several hundred Hz, 
so there seems to be something limiting the accretion- 
induced spin-up. A limiting mechanism that has been 
suggested is the emission of gravitational waves Danin], 
which would result in steady-state GW emission where 
the accretion torque is balanced by the radiated angu¬ 
lar momentum. For a discussion of alternative expla¬ 
nations, see n IS]. The resulting torque-balance GW 
amplitude increases with the observed X-ray flux (inde¬ 
pendent of the distance of the system, e.g. see Eq.(4) in 
m), therefore the X-ray brightest LMXB Scorpius X-1 is 
often considered the most promising GW source within 
this category mini- 

several searches for GW signals from Scorpius X-1 have 
been performed (without any detections) on data from 
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initial LIGO [TTI - H3] , and several new pipelines have been 
developed and are currently being tested in a Scorpius X- 
1 Mock Data Ghallenge (MDG) [lOj . So far, all current 
methods fall into either one of two extreme cases: highly 
coherent with a short total time baseline (6 hours in [TT] , 
10 days in USD, or highly incoherent with a long total 
baseline but very short (^hours) coherent segments [121 

nmn]. 

In order to increase sensitivity beyond these methods, 
and to be able to effectively absorb large amounts of com¬ 
puting power (such as those provided by Einstein@Home 
[HI El), it is necessary to extend the search approach 
into the realm of general long-segment semi-coherent 
methods, by stacking coherent J^-statistic segments nn- 
[30] . Such methods have already been employed over the 
past years for both directed and all-sky searches for GWs 
from isolated NSs iniEIlES], but they have not yet been 
extended to the search for GWs from binary systems. 

Key ingredients required for building a semi-coherent 
“StackSlide” search are the coherent and semi-coherent 
parameter-space metrics na [33H3S]. The coherent bi¬ 
nary GW metric was first analyzed in [26], ^^nd this was 
further developed and extended to the semi-coherent case 
in m- Here we will largely follow the approach of EH, 
considering only low-eccentricity orbits, and focusing on 
two different regimes of either long coherent segments, or 
short segments compared to the orbital period, both of 
which admit simple analytic results. 


Plan of this paper 


Sec. |TT| provides a general introduction of the con¬ 
cepts and notation of semi-coherent StackSlide methods, 
parameter-space metric and template banks. Secs. IIT Vj 
build on the work of [27], rectifying some of the results 
and extending them to the case of a general orbital refer- 





2 


ence time. In Sec. VI we present the results of extensive 
Monte-Carlo software-injection tests comparing the met¬ 
ric predictions against the measured loss in statistics due 
to parameter-space offsets. Sec. |VII| applies the theo¬ 
retical template-bank counts to compute optimal Stack- 
Slide setups for a directed search for Scorpius X-1 and 
estimates the resulting achievable sensitivity. Finally, 
Sec. |VIlI| presents concluding remarks. 

Notation. Throughout the paper we denote a quantity 


Q as Q when referring to the coherent case, and as Q 
when referring to the semi-coherent case. 


II. BACKGROUND 

The strain h in a given detector due to a continuous 
gravitational wave is a scalar function h(t;A,X), where 
t is the time at the detector. The set of signal parame¬ 
ters A denotes the four amplitude parameters, namely the 
overall amplitude h^, polarization angles cos l and ^/>, and 
the initial phase (jjo- The set of phase-evolution parame¬ 
ters A consist of all the remaining signal parameters af¬ 
fecting the time-evolution of the CW phase at the detec¬ 
tor, notably the signal frequency /, its frequency deriva¬ 
tives (also known as “spindown” terms), sky-position of 
the source, and binary orbital parameters. We will look 
in more detail at the phase model in Sec. [ml 

The observed strain x{t) in the detector affected by 
additive Gaussian noise n{t) can be written as x{t) = 
n(t) h{t;A,X), and the detection problem consists in 
distinguishing the (pure) “noise hypothesis” of h-o = 0 
from the “signal hypothesis” hg > 0. 


by the required computing cost. It was shown that semi- 
coherent statistics generally result in better sensitivity at 
equal computing cost (e.g. see mm)- 


B. Semi-coherent detection statistic 


Here we focus on one particular semi-coherent ap¬ 
proach, sometimes referred to as StackSlide, which con¬ 
sists of dividing the total amount of data Tobs into N 
segments of duration AT, such that Tobs = N AT (in 
the ideal case of gapless data). The coherent statistic 
J-i{x-, A) is then computed over all segments £ = 1... N, 
and combined incoherently by summing 

N 

.F(a.;A) = ^.T,(x;A). (2) 

e=i 

This statistic follows a (non-central) y^-distribution with 
47V degrees of freedom and a non-centrality parameter 
p^{A, As; A), which is found to be given by 

N 

= (3) 

e=i 


in terms of the per-segment coherent SNRs pg. Note that 
(contrary to the coherent case) 'p can not be regarded as 
a semi-coherent SNR. The expectation for T is 


E 


2T(x-,X) 


— 47V -|- p^{A, As; A). 


(4) 


C. Template banks and metric mismatch 


A. Coherent detection statistic 

As shown in [55], matched-filtering the data x{t) 
against a template h{t-,A,X), and analytically maximiz¬ 
ing over the unknown amplitude parameters A, results 
in the coherent statistic T'(a;;A), which only depends on 
the data x and on the template phase-evolution param¬ 
eters A. This statistic follows a (non-central) distri¬ 
bution with 4 degrees of freedom, and a non-centrality 
parameter ^(M, As;A), which depends on the signal pa¬ 
rameters {A, As} and the template parameters A. The 
quantity p (which depends linearly on Tiq) is generally 
referred to as the signal-to-noise ratio (SNR) of the co¬ 
herent T-statistic. The expectation value of this statistic 
over noise realizations is given by |55| 


In order to systematically search a parameter space 
V using a statistic J-'(a;;A) (which here can refer either 
to the coherent T or the semi-coherent T), we need to 
select a finite sampling {At} C T of the parameter space, 
commonly referred to as a template bank, and compute 
the statistic over this set of templates, i.e. {T(a;; At)}. 

A signal with parameters As G T will generally not fall 
on an exact template, and we therefore need to charac¬ 
terize the loss of detection statistic T as a function of the 
offset (5A = A — As from a signal. This is generally quan- 
tihed using the expected statistic E [2T] of Eqs. ([T) and 
Q, namely (by removing the bias 47V) as the relative 
loss in the non-centrality 


_ p'^{A, As; As) - p'^{A, As; A) 


(5) 


E 


2T{x-,X) 


4 -f {A, As; A). 


( 1 ) 


In the case of unknown signal parameters As S T within 
some parameter space V, the number of required tem¬ 
plate At that need to be searched is typically impracti- 
cally large. The maximal achievable sensitivity by this 
coherent statistic ^( 0 ;; A) is therefore severely limited [55] 


^ This form of the “ideal” StackSlide statistic is not directly used 
for actual searches. For reasons of computational cost, the coher¬ 
ent statistics across segments are computed on coarser template 
banks, and are combined on a fine template bank by summing 
across segments using interpolation on the coarse grids. This 
is discussed in detail in m> but is not relevant for the present 
investigation. 
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which defines the measured J^-statistic mismatch func¬ 
tion. This mismatch has a global minimum of /tq = 0 at 
A = As and is bounded within /tq € [0, !]■ 

We now use two standard approximations: (i) Taylor- 
expand this up to second order in small offsets (5A, and 
(ii) neglect the dependence on the amplitude parameters 
A (the effect of which was analyzed in detail in [29)'). 
which leads to the well-known phase metric pij, namely 


place templates in order to limit the maximal (relative) 
loss in detection statistic, which can be written as 

max min/r(As; At) < Mmax ■ (H) 

As At 

This states that the worst-case mismatch to the clos¬ 
est template over the whole template bank should be 
bounded by a maximal value /imax- Note that each tem¬ 
plate At “covers” a parameter-space volume given by 


/ro(Al, As; A) « ^(As; A) , where (6) 

/r(As;A) =5 ,,(As)(5AW, (7) 

with gij a positive-definite symmetric matrix, and im¬ 
plicit summation over repeated indices *, j is used. The 
quality of this approximation will be quantified in a 


Monte-Carlo study in Sec. VI The metric mismatch 
pL has a global minimum of /x = 0 at A = As and 
(contrary to /xq) is semi-unbounded, i.e. /x S [0,oo). 
Previous studies testing p, against /xq PM5T] for all¬ 
sky searches have shown that the phase-metric approx¬ 
imation Eq. §) works reasonably well for observation 
times Tobs C>(lday) and for mismatch values up to 
/X < 0.3 — 0.5. Above this mismatch regime, the metric 
mismatch /x increasingly over-estimates the actual loss 
/Xq. For example, in |29| an empirical fit for this behavior 
was given as /xo(/x) « /x — 0.38/x^. 

The metric formalism was first introduced in [23l [24] 
(in the context of compact binary coalescence signals), 
and was later applied to the CW problem in [^, where 
it was shown that the coherent phase metric can be 
expressed explicitly as 


gy(A) = dj(j){X)) - (9*())(A)) {dj(l){X)) , (8) 

where dtcj) = d(j)/dX'‘, with (j){i\X) the CW signal phase, 
and (...) denotes the time average over the coherence 
time r, i.e. {Q) = Qi^) 

The corresponding semi-coherent phase metric gij was 
first studied in CHI and was found to be expressible as 
the average over all per-segment coherent phase metrics]^ 


9ij{^) — 9ij{^) ) ( 9 ) 

where we defined the average operator over segments as 
_ 1 ^ 

The metric is useful for constructing template banks, by 
providing a simple criterion for how “close” we need to 


/x(At; A) = gij[Xt) <5A* 5X^ < /x^ax , (12) 

which describes an n-dimensional ellipsoid, where n is 
the number of parameter-space dimensions. The tem¬ 
plate bank can therefore be thought of as a covering of 
the whole parameter space with such template ellipsoids, 
such that no region of V remains uncovered. One can 
show [25l |32l [33] that the resulting number of templates 
N in such a template bank is expressible as 


W = 0„/x-:/2 [ ydrt^d"A, (13) 

Jv 

where n is the number of template-bank dimensions, and 
9n is the center density (also known as normalized thick¬ 
ness) of the covering lattice, which quantifies the number 
of lattice points per (metric) volume for a unit mismatch. 
The center density is a geometric property of the lattice, 
and for the typical covering lattices used here is given 
by [32] 


f 2-"n”/2 


0. = 

n(n-|-2) 

12(n-|-l) 

Vn -1- 1 


1 "/2 


for , 
for A! . 


(14) 


An important caveat for using Eq. (13) is that the met¬ 


ric determinant and integration must only include fully 
resolved template-bank dimensions, which require more 
than one template to cover the extent AA* of the param¬ 
eter space along that dimension. We can estimate the 
per-dimension template extents at given maximal mis¬ 
match ^niax I 

(5V = 2 V/imax [<?-']" , (15) 


where \g~^X^ are the elements of the inverse matrix to the 
metric gij. We can therefore define the metric template- 
bank “thickness” along dimension i in terms of the cor¬ 
responding effective number of templates along that di¬ 
rection as 





(16) 


For dimensions A* with J\f\i < 1 we must exclude 
this coordinate from the bulk template-counting formula 
Eq. ( [T3| ), as it would effectively contribute “fractional 
templates” and thereby incorrectly reduce the total tem¬ 
plate count. 


^ This assumes constant signal SNR As; As) over all segments 

£, which should be a good approximation for stationary noise 
and sufficiently long segments, such that diurnal antenna-pattern 
variations have averaged out. 


^ Note that the extra factor of 2 here compared to [7] is required 
to account for the total extent of the ellipse, not just the extent 
from the center. This can also be seen by considering that the 
center density = 1/2 in Eq. ( |13[ l for the 1-dimensional case. 
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III. THE BINARY CW PHASE 


A. The general phase model 


The general CW phase model assumes a slowly 
spinning-down NS with rotation rate and a 

quadrupolar deformation resulting in the emission of 
CWs. The phase evolution can therefore be expressed 
as a Taylor series in the NS source frame as 

= 27r [/ (r - tref) + ^/(t - Uetf + ..■], (17) 

where tref denotes the reference time and /,/,/,■•■ are 
the CW frequency and spindown parameters. These are 
defined as 






r = tret 


(18) 


where cr is a model-dependent constant relating the in¬ 
stantaneous CW frequency /(r) to the NS spin fre¬ 
quency v{t). For example, for a quadrupolar deforma¬ 
tion rotating rigidly with the NS (“mountain”) we have 
tj = 2 [31 ini I31H3S], while for r-modes a « 4/3 jH[371 |3H] 
and for precession a » 1 [311110]. 

In order to relate the CW phase in the source frame 
to the phase />(tarr) in the detector frame needed for 
Eq. Q, we need to relate the wavefront detector arrival 
time tarr to its source emission time r, i.e. T(tarr)) such 
that ^!)(tarr) = />’’‘^'^(;l(tarr))- Neglecting relativistic wave- 
propagation effect^ (such as Einstein and Shapiro de¬ 
lays, e.g. see mi if for more details), we can write this 
as 


^(tarr) — tarr 4“ 


7'(tarr) ’ n d R{t) 


(19) 


where r is the vector from solar-system barycenter (SSB) 
to the detector, c is the speed of light, d is the (gener¬ 
ally unknown) distance between the SSB and the binary 
barycenter (BB), R is the radial distance of the CW- 
emitting NS from the BB along the line of sight, where 
i? > 0 means the NS is further away from us than the 
BB, and n is the unit vector pointing from the SSB to the 
source. In standard equatorial coordinates with right as¬ 
cension a and declination d, the components of the unit 
vector n are given by n = (cos a cos <5, sin a cos <5, sin 5). 

The projected radial distance R along the line of sight 
can be expressed (e.g. see [43]) as 


R = y sini + v), (20) 


where y is the distance of the NS from the BB, i is the in¬ 
clination angle between the orbital plane and the sky, uj is 


the argument of periapse, and v is the true anomaly (i.e. 
the angle from the periapse to the current NS location 
around the BB). We further approximate the orbital mo¬ 
tion as a pure Keplerian ellipse, which can be described 
as 


y{v) 


a (1 — e^) 

1 -|- e cos V ’ 


( 21 ) 


in terms of the semi-major axis a and the eccentricity e. 
The ellipse can be written equivalently in terms of the 
eccentric anomaly E, namely 


y{E) = a(l — e cosE), (22) 


and the dynamics is described by Kepler’s equation, 
namely 


T -tp = — {E - e sin E) , 

ZTT 


(23) 


which provides a (transcendental) relation for E{t). 
Combining Eqs. (20), (21) and (22), we can rewrite the 
projected radial distance R in terms of E, namely 


R 

— 

c 


sinw(cosE — e) -I- cosa;sinE\/l — , (24) 


where we defined Op = a sini/c. Combining this with 
Eq. (23) fully determines (albeit only implicitly) the func¬ 
tional relation R{t) required for the timing model of 

Eq. ([^. 

Dropping the unknown distance d to the BB (which is 
equivalent to re-defining the intrinsic spindown parame¬ 
ters), and defining the SSB wavefront arrival time tsSB 
as 


IsSB(larri^) — ^arr 4” 


r{ta 


we can rewrite the timing relation Eq. (191 as 

, , R{t) 

T{tsSB) = tsSB -• 


(25) 


(26) 


Here we are interested only in binary systems with 
known sky-position ft, and we can therefore change in¬ 
tegration variables from t to tsSB in the phase for the 
metric integration of Eq. (|^. Furthermore, given that 
the Earth’s Rpmer delay is bounded by |r • n/c\ < 500 s, 
and |dtsSB/dt| ^ 10“^, this difference will be negligible 
for the metric and can be dropped. This is equivalent to 
effectively placing us into the SSB, which is always pos¬ 
sible for known n. In order to simplify the notation, we 
now simply write t = tssB ■ 

PluggijM this timing model into the phase of Eq. (0, 
we obtaiio 


(f>{t) « 27r 


/ At 






(27) 


^ These effects are taken into account in the actual matched- 

filtering search codes, but are negligible for the calculation of ® The sign on it (and equivalently on Up) here agrees with Eq.(2.26) 
the metric. in but differs from Eq.(2) in m which is incorrect. 
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with At = t — tref- The binary systems we are interested 
in have semi-major axis Op of order 0(1 — 10 ) s, and bi¬ 
nary periods P of order of several hours. Hence, the 
change in E, and therefore R{E), during the time R/c 
will be negligible, and so we can approximate E{t) « 
E{t), namely 


P 

t — tpKi — {E — e sin E) . (28) 

27r 

For the purpose of calculating the metric using Eq. Q, 
we can further approximate the phase model in the fol¬ 
lowing standard way (e.g. see ESEO]) as 


(j){t) ^ 27T 


{fAt+^fAt^+ ...)- 


(29) 


where we expanded the factors (At — R/c)^ and kept 
only the leading-order terms in At, keeping in mind that 
T > O(days) and R/c ~ O(seconds), where T de¬ 
notes the coherence time. In the last term we replaced 
the instantaneous intrinsic CW frequency as a function 
of t by a constant parameter /, namely we approximated 
/«(/-!- /At -I- ...). This “frequency scale” / of the 
signal could be chosen as the average (or the largest) 
intrinsic CW frequency of this phase model over the co¬ 
herence time T. Given that this only enters as a scale 
parameter in the metric, and the changes of intrinsic fre¬ 
quency over the observation time will generally be small, 
this introduces only a negligible difference. 


B. The small-eccentricity approximation (e <C 1) 


We follow the approach of m and consider only the 
small-eccentricity limit of the metric. Hence, by Taylor 
expanding Eqs. (241 and ( p^ up to leading order in e, 
i.e. inserting E{t) = Eo{t) -I- eEi{t) . into Kepler’s 
equation Eq. (28), we obtain 


Eo{t) = n{t - tp), (30) 

Ei(t) = sin£'o(t), (31) 


where H ^ is the mean orbital angular velocity. Plug¬ 
ging this into Eq. (24), we obtain the Rpmer delay of the 
binary to leading order in e as 


R 

c 


sin ip{t) -\- 


K 

2 


sin 2ip{t) 


^cos2tjj{t) , 


(32) 


where a constant term —Zaprj/2 was omitted, which is 
irrelevant for the metric. We use the standard Laplace- 
Lagrange parameters defined as 

K = e cos(w), (33) 

?7 = e sin(a;), (34) 

and the mean orbital phase 

-tjj^t) = H (t - tasc) , (35) 


measured from the time of ascending node tasc, which 
(for small e) is related to tp by [4T] 

tasc = ^ ’ (3®) 

and which (contrary to the time of periapse tp), remains 
well-defined even in the limit of circular orbits. 

The small-eccentricity phase model is there¬ 
fore parametrized by the 5 binary parameters 
{op, tasc) ??} (referred to as the “ELLl” model in 
mHHDj and in the circular case (e = 0 ) this reduces to 
the 3 binary parameters {op, taso with k = 77 = 0. 


IV. THE BINARY CW METRIC 


A. Phase derivatives 


Following [Ml \22\ we restrict our investigation to 
(approximately) constant-frequency CW signals, i.e. 
y(fc>o) _ This is motivated by the assumed steady- 
state torque-balance situation in LMXBs, which are our 
main target of interest. However, the corresponding fluc¬ 
tuations in the accretion rate are expected to cause some 
stochastic frequency drift, and one will therefore need 
to be careful to restrict the maximal coherence time in 
order to limit the frequency resolution. This will be dis¬ 
cussed in more detail in the application to Scorpius X-1 


in Sec. VH The total phase-evolution parameter space 


considered here is therefore spanned by the following 6 
coordinates 


A = {/, Op, tasc, K, 77 } . 


(37) 


The small-eccentricity phase model used here can now be 
explicitly written as 

/ At -/op [sin/;-b ^sin2V) - ^ cos 2 V 7 j , 

ZTT L Z Z J 

(38) 

with if = ifiyt) given by Eq. (35). The frequency param¬ 
eter / is treated as a constant in the phase derivative 
df(j) (which corresponds to neglecting the small correc¬ 
tion Op <C At ~ T), but numerically we have / = / in 
the present constant-frequency case without spindowns. 
Hence, we obtain the phase derivatives 


27T 

27T 

2tt 


^ At 
= -/ 


sin if -i- sin 2if — ^ cos 2if 


= /opH [cost/-I-kcos 27/)-I-?7sin27/)] , (39) 


d (j) 

—— = -fap{t - ta.sc)[cos if K cos 2if -I- 77 sin 2 V!], 

ZTT 

d^(f 1 J 

= - -fapsm2if, 

^fapCOs2if, 
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which are inserted into Eq. ([^ in order to obtain the 
coherent metric. The semi-coherent metric is obtained by 
averaging the coherent metrics over segments according 
to Eq. ([^. 

The resulting analytic expressions for these metrics in 
the general case are quite uninstructive and unwieldy, 
while it is straightforward to compute them numerically 
for any case of interest. However, as noticed in previous 
investigations (23 HZ], it is instructive to focus on two 
limiting regimes that yield particularly simple analytical 
results, namely the long-segment limit (LS) where AT 3> 
P, and the short-segment limit (SS) where AT <?; P. 

When taking these limits on the metric gij , in order to 
decide whether a particular off-diagonal term is negligible 
or not, it is useful to consider the diagonal-rescaled metric 


DN _ 9ij 
"ij — 


\/9ii 9jj 


(40) 


which is dimensionless and has unit diagonal, = 1. 
In this rescaled metric we can then naturally neglect off- 
diagonal terms if they satisfy ^ 1, as their corre¬ 
sponding contribution to the metric mismatch of Eq. (12) 
will then be negligible. 


B. The long-segment (LS) regime (AT » P) 


where the only nonzero components of the Jacobian of 


this coordinate transformation Eq. (41) are 


d\^ _ 

'W ~ ’ 


and 


dt: 


an 


= -p 


27r 

■ 


(44) 


In the long-segment limit, these discrete steps P = 2tt/Q, 
are assumed to be small compared to the segment length 
AT. We can therefore choose the time of ascending node 
to be approximately at the segment midtime tmid, he. 
we consider the special gauge choice « tmid, in which 
the phase derivatives Eq. (391 are time-symmetric around 


tmid, which simplifies the metric calculation. Introducing 
the time offset 


Ania — i'mid ^asc , (^^) 

we therefor consider first the special case = 

5l)®(Ama = 0). 

Keeping only leading-order terms in H AT 3> 27r, and 
up to first order in e (i.e. first order in k and rj), we obtain 
the diagonal metric 

ffmid.ii = X (46) 

' AT^ -2 (/«pAT)^ ifa^ {Ja^r 

g ; / I U p 1 ; ^2 ’ 4 ’ 4 ’ 


1. Coherent metric 


As noted in [27|, it is convenient to use the discrete 
gauge freedom in the choice of the orbital reference epoch 
tasc, as we can choose the time of ascending node during 
any orbit. Therefore we can add any integer multiple p 
of a period and redefine 

QiTT 

4c=^asc+P-^, PeZ, (41) 


without changing the system. An (infinitesimal) offset 
(or “uncertainty”) Jtasc changes under such a transfor¬ 
mation in the presence of an (infinitesimal) offset on the 
period, i.e. if Sfl 0, namely 

, 27r 

<^iasc = ^^asc “ Pt^ dCl . (42) 


All other coordinate offsets are unaffected by this change 
of orbital reference epoch. This needs to be carefully 
taken into account in the metric when redefining tasc- 
Namely, given that this is a pure “relabeling” of the same 
physical situation, the corresponding metric mismatch 
must be invariant, i.e. 


p=6X^g'^,SX 


'b 


=SX^ 


ax 


ax 


% 


ax 

=5X g,j SX^ , 


9mb 


axi 


6X^ 


(43) 


in agreement with the result Eq. (30) of |?7j . 

The general case of A^a 0 can be obtained by apply¬ 
ing the Jacobiaij^of Eq. (44) with A^a = p'^.-nj^-, and so 
^'asc = ^asc + Ama and at^j aXt = -Ama/fi, resulting in 
the general coherent long-segment metric (7]^®(Ania) with 
components 


TtLS _ 2 

9ff - ^ 


AT2 


9t. 


yapO 

'LS 


3 ’ 

= P , 

= (/opH)^ , 


9nn — 271" (/op) 

5nLc = = -27r^ {fapff n A„ 


AT2 


9^^ 


TT 


PJ \ ^2 

-ZV.2 

2 


ma 


= 9vv=Y ’ 


(47) 


while all other off-diagonal terms are (approximately) 
zero in this limit. The nonzero off-diagonal term shows 
that in general (i.e. for A^a 0) there are correlations 
between offsets in tasc and H, introduced by the relation 
Eq. (42). This generalizes the result in [27] to arbitrary 


choices of orbital reference epoch t^ 


A direct but somewhat more complicated calculation of the met¬ 
ric at Ama 7 ^ 0 yields the same result 
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Z. Semi-coherent metric 


C. The short-segment (SS) regime (AT <P) 


In order to compute the semi-coherent metric, we only 
need to average the per-segment coherent metrics of the 
previous sections over segments according to Eq. 

Here it is crucial to realize that all segment metrics must 
use the same coordinates A in order for this averaging ex¬ 
pression to apply, in particular we can only fix the gauge 
on tasc once for all segments, and so generally Ama,^ ^ 0 . 
This means we cannot use the special form of of 
Eq. (46 1 for the average over all segments, as done in 


which would result in the erroneous conclusion that 
the semi-coherent metric would be identical to the 
coherent one. Instead, using Eq. ([^, we find by averag¬ 
ing the per-segment metrics 5 ^® = g^®(Ai„a,£): 


^LS _ 2 

9ff - ^ 


AT 2 


^LS 

JdpQp 


-LS 


3 ’ 

= P , 

= (/apH )2 , 


(48) 


dan = 27r^ {fopf -F A^a^ 


-LS 


= dtln = -27r^ ifopf n Ar 


^LS 

ciKK 


= 9l;l = \{fa^r 


It is interesting to consider the ideal case of N regularly- 
spaced segments without gaps, i.e. Tobs = N AT, and 


Ania,r (^mid f'asc) 4“ 

resulting in the first two moments 


N+1 


AT, 


Ania — ^mid ^a 


A^a = A„.a +-(iV^-l)AT^ 


(49) 

(50) 

(51) 


where tmid is the average mid-time over all segments. In 
this case we can write the variance of {Ama,^} as 

I 


var [A„,a] = A^a - A„,a =-(7V^-l)AT^ (52) 


If we use the g auge freedom and choose tasc ~ ^mid; then 
Ama ~ 0, and A^a ~ — 1) AT^/12, so that in this 

case the metric is again diagonal, and the only component 
different from the coherent metric is 

ghn = (A AT)2 cx t2,, , (53) 

corresponding to a refinement in the fl coordinate. 
Therefore the semi-coherent template-bank spacing in H 
needs to be finer by a factor of N compared to coherent 
spacing in order to achieve the same mismatch. The ex¬ 
istence of this metric refinement in 17 in the long-segment 
limit had already been noticed earlier [45], and was also 
used in the discovery of a binary pulsar in Fermi LAT 
data US). 


1. Coherent metric gt®® 


In the case of very short coherent segments compared 
to the period, i.e. flAT 27r, as pointed out in m, 
strong parameter-space degeneracies render the metric 
in these coordinates quite impractical for template-bank 
generation, and we therefore follow their approach of 
Taylor-expanding the phase around the mid-point tmid 
of the coherent observation time, namely 

n 

(j){t) = 27r ^ ^ (t - tmid)'", (54) 

fe=i 


omitting a constant term. 

The M-coordinates {uk} are defined as the k-th time- 
derivatives of the phase at tmid; namely 




(55) 


Note that the phase model of Eq. (54) is formally identi¬ 


cal to the phase in terms of the frequency and spindowns 
of an isolated NS (with tmid playing the role of 
the reference time tref), as can be seen in Eq. (29). The 
corresponding “spindown metric” has a well-known ana¬ 
lytical form, which can be expressed most conveniently 
in terms of the rescaled dimensionless ri-coordinates. 


Vk = 2n 


Uk 

k\ 



(56) 


resulting in the z;-coordinate metric 




{k -\- k' — I) mod 2 
k -\- k' -\-l 

[{k — 1) mod 2] [{k' — I) mod 2] 


n 

0 

1 

5 

0 

...\ 

0 

4 

0 

8 


45 

105 


1 

5 

0 

1 

7 

0 


0 

8 

0 

16 


105 

225 


V 
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(57) 


correcting the incorrect expression Eq. (21) of [ 27 ] . 

Applying Eq. (55) to the small-eccentricity phase of 
Eq. (38), we obtain the following expressions 


Ui = —fOpfl [cos Ipjn K cos -|- T] sin 2'0m] + /; 

U 2 = fapCl^ [sin + 2k sin 2il}^ — 2r] cos 2 '!/;m], 

U3 = [cos'i/'m + 4 KC 0 s 2 '!/)m + 4 ? 7 Sin 2 ' 0 m], (58) 

Ui = —fopCl'^ [sin^m -I- Sk sin 2'(/)i„ — 877 cos 2'0m], 

U5 = — /apn^[cOS V'm + 16KC0s2'!/;m + 1677 sin 2'0ni], 
uq = fapCl^ [sin -|- 32k sin 2tp^ — 32g cos 2ip ^], 




























with = '0(tmid), which differ in the sign of Op com¬ 
pared to Eqs. (A2-A5) of [27], due to the sign error on R 
in the phase model discussed in Sec. |III A| 

Note that for the general elliptical case we expect 6 in¬ 
dependent u-coordinates m(A) and 4 in the circular case. 
The nonlinear coordinate transformation Eqs. (58) can 


be analytically inverted to obtain A (it), as shown explic¬ 
itly in Appendix |A| 

We follow in estimating the range of validity of this 
Taylor approximation up to order n by considering the 
order of magnitude of the first neglected term of order 
n -|- 1, namely 


|A^| ~ 27r 


'^n +1 


27r fa 


AT A 

(n-hl)! / 

- (7ra:)”+i 


n+1 


(n + 1)! 


(59) 


where we defined the fraction x of an orbit sweeped dur¬ 
ing the segment duration, i.e. x = nAT/(27r) = AT/P. 
In order to link this phase error to a mismatch, we ob¬ 
serve that here A(j) ~ Ai;„+i of Eq. (56), and the corre¬ 


sponding metric element of Eq. (57) gives (n+i) 


271-|-3 


for n ^ 1 and therefore 

_Af^_ - ,2 (7rx)^("+^) 

^ 2n -h 3 (2n -h 3) (n -H 1)!^ 


(60) 


Plugging in typically values used in the numerical tests 
later, e.g. / ~ 500 Hz, Op ^ 3 s, n = 6, we see that 
the mismatch grows very rapidly from /i ^ 1 at a; < 
0.25 to /I ~ 1 at a; ~ 0.35. We therefore expect this 
approximation with n = 6 to typically hold for segment 
durations up to a 1/4 to 1/3 of an orbit. 


2. Semi-coherent metric 


While the segments AT are assumed to be short com¬ 
pared to the period P, we only consider the case where 
the total observation time is long compared to P i.e. in 
addition to AT <C P we also assume Tobs = NAT ^ P, 
which also implies N ^ 1. 

We cannot directly use the coherent per-segment met¬ 
rics in li-coordinates to compute the semi-coherent met¬ 
ric as an average using Eq. ([^, because a fixed phys¬ 
ical parameter-space point A would have different u- 
coordinates in each segment. 

Instead, as first shown in 113, one can go back to phys¬ 
ical coordinates and use the fact that A ^ 1 to replace 
the discrete sum in Eq. 0 over segments by an integral 
over time, namely 


N 


5^® = 

Af»l 1 


r^mid“t-T)bs/2 


T 


obs .7t„id-Tobs/2 




(61) 


in terms of coherent per-segment metrics 5^®(tmid,£) ex¬ 
pressed as a function of the segment mid-times tmid.i, 
and where t^id is the mid-time of the whole observation. 
In order to compute this, we first calculate the coherent 
metric of Eq. ^ for each segment —AT/2, tmid,£ + 
AT/2] by using the phase derivatives of Eq. (39). As be¬ 
fore, we keep only first-order terms in e (i.e. first-order in 
K, ly), then Taylor-expand the results in HAT ^ 27r up 
to second order, obtaining ^^(tmid.t)- This is integrated 
over the total observation time according to Eq. (61), 


and only leading-order terms in HTobs ^ 27r are kept. 
We only keep off-diagonal terms where the correspond¬ 
ing diagonal-rescaled elements Eq. (40) are not ^ 1. The 


resulting semi-coherent short-segment long-observation 
limit metric elements are found as 


9ff 


AT^ 


JQpClp 


= -{nATfp, 


y^ascia 


= —{^lATf 


=-(HAT) 2 (/ap) 


(62) 


rp2 

2 I ^obs 
12 




= Cta = -y(f^AT)^(/ap)^HA„ 


JK.K 


= ^“ = ^(HAT)^(/ap)^ 


which generalizes the result in HZ] (for gauge choice 
^asc ~ ^mid , i.e. A,„ia ~ 6) fo the general case of A^a 7 ^ 0- 
Note that a more general form encompassing the semi- 
coherent metric in both the short-segment and long- 
segment limits has recently found in m- 


V. NUMBER OF TEMPLATES 


In order to express the explicit template counts based 
on Eq. (|13|), we will assume a simple parameter space V 


bounded by / S [/min;/max]) Op S [®Pmin’®Pmax] ’ 
[^ascmin; 


^-ascmm; ^'ascmaxj; 
^ ^ [o^min;Oimax 




e e 


H G 

. Using Eq. (R^, (pMj) we see that 


[Omin;emax]; and 


dndr] = 


ededu). 


(63) 


Observing that g^K = g-riv^ further obtain /i = 

-\- grjn drf = gK.K {de^ -I- duP'), which implies 

9ee ~ ShlK ; guiUJ — e gKK. ; geu) — 0 • (^4) 

Eor convenience of notation, for the following template- 
count expressions, we define the shorthand 


\Q] = Q 


-Qn 


(65) 


max 
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A. Templates in the long-segment regime 

The determinant of the semi-coherent metric of 
Eq. (48) only differs from that of the coherent metric 
of Eq. (471 by the presence of a rehnement factor 7 , i.e. 


det = 7 ^ det g 


:LS 


with 


7 = Wl-tl 2 


var [A 


AT 2 


maj gaptess ^ 


( 66 ) 

(67) 


where the last equality holds in the special case of N 
segments without gaps, as seen from Eq. (52). This re¬ 


finement only affects the f2-dimension. In particular, for 
the gauge choice A^a = 0 , we see 


-^LS 2 -^LS 

9Tnid,flfl ^ ' 


( 68 ) 


The determinant Eq. (661 is seen to be independent of the 
gauge choice on tasc (be. on A^a), and we can therefore 
write the template volume density in factored form as 


y/det g^^ dTX = 




dy 


(69) 


i=l 


This independence of A^a in the volume element may 
seem surprising at first, as one would expect |48| the 
parameter-space uncertainty on t^sc to grow with in¬ 
creasing separation from the observation time tmid = 
tasc + Ama, if there is any uncertainty d^l in the mean 
angular velocity, as also discussed in Sec. |IVB| However, 
this is only the case if the uncertainty in is not re¬ 
solved by the template bank, when one observes indeed 
a stretching of the effective (projected) uncertainty in 
tascj resulting in an increase in required templates 
On the other hand, if the O —tasc parameter space is fully 
resolved by the metric [as assumed in Eq. ([ 6 ^], there is 
no increase of templates. Changing Ama only deforms 
the O — tasc parameter space in a volume-preserving way 
by shear according to Eq. (411. Given that any offset in 
tasc and O results in a deterministically shifted offset 
there is no loss of information, and therefore no increase 
in the number of templates. Alternatively, one can con¬ 
sider the orbital epoch t^sc with its uncertainty hxed at 
its original measurement epoch, and account for the off¬ 
set Ajna in the metric as done in Sec. HYl This results 
in changing metric correlations between tasc and while 
leaving the ellipse volume unchanged. 

The expression for the number of templates per dimen¬ 
sion of Eqs. (IT^ and (flbl) simplifies to 


= \ 




(70) 


with individual components in coordinates A* = 
{/, Op, tasc, e} given by 




, 7 Ar /flp 

v^,/’/“P ,^“P v^, 2 


( 71 ) 


noting that gee = = dw gu,u, = gee- In prac¬ 

tice, one often encounters cases where the parameter un¬ 
certainty along some of these dimensions is smaller than 
the metric resolution, such that a single template covers 
the whole extent of the parameter space along that di¬ 
rection. In this case the corresponding coordinate contri¬ 
bution to the template density would result in fractional 
templates (and therefore underestimating the number of 
templates) and must not be included in Eq. (13), as dis¬ 
cussed in Sec. Eel 

The number of templates over the full 6 D parame¬ 


ter space {/, Op, tasc) e, w} according to Eq. (13) is ob¬ 
tained as 




’ 7 AT' 


Mmax 


36 OV 2 


^ \r 


\ta 




(72) 

In the application to Scorpius X-1, a few special cases will 
be of interest, namely when one or more of the uncertain¬ 
ties in orbital parameters are smaller than the template 
extent. 

For sufficiently well-estimated orbital angular velocity 
n = rio (such that Nn < 1 ), the template count for the 
corresponding 5D template bank over {/, Op, taso e, w} is 
found as 


= 


■AT 


5/2 

/^max 


40^3 


IIo \f] I Op] |tasc] |e^] |w] , (73) 


while in the 4D case of well-estimated e and uj (or equiv¬ 
alently, for circular orbits) {/, Op, tasc, bl}, we obtain 




O4 qAT^ 
36 


Mmax 


1/"^] |«p] |tasc] 


(74) 


and hnally in the 3D circular case with well-determined 
D = Dg, e and uj we hnd the template count over the 
remaining parameter space {/, Upjtasc} as 


= 


6»3 tt^AT 

A 


flo \P 


|ta 


(75) 


We note that Eqs. (ln])-([^ are valid for the semi- 
coherent case with general refinement 7 , as well as for 
the coherent case with 7 = 1 . 


B. Templates in the short-segment regime 


Due to the nonlinear transformation Eqs. (58) from 
physical parameters A into u-coordinates, the physical 
parameter space V would be described by complicated 
integration boundaries in u. In addition, the number of 
dimensions in u-coordinates that need to be included in 
the template bank is a non-trivial function of the signal 
parameters and the mismatch /imax (e.g. as seen later 
in Fig. [^. These effects would dominate the expression 
for the number of templates Eq. (13), and it is therefore 


not clear whether a useful closed analytical form can be 
given. 
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Using the metric volume density for the semi-coherent 
short-segment metric of Eq. (62), we obtain the fol¬ 
lowing number of templates Eq. (131 for the full 6D tem¬ 
plate bank 


= 




i rj. 


obs 


Mmax 


90720V6 

while in the 4D circular case we hnd 


f] I Op] |tasc] 1^^^] |e^] Iw] 


= 


04 ^^AT^Tobs ,X4 


Mmax 


2160^6 


n 


l^asc] 




(76) 

(77) 


VI. NUMERICAL TESTS OF THE METRICS 

In this section we present numerical tests performed 
on the parameter-space metrics derived in the previous 
sections. These tests consist in comparing the predicted 
metric mismatches fi of Eq. 0. against measured J"- 
statistic mismatches no of Eq. ([^ obtained via signal 
software injections. The quantity used for this compari¬ 
son is the (symmetric) relative difference [2M^ , defined 
as 


Choice of signal parameters {A, As} 


The random signal amplitude parameters A are chosen 
as follows: the scalar amplitude /ig plays no role for the 
metric and is fixed to hg = 1. The inclination angle is 
drawn from a uniform distribution in cosi S [—1,1], the 
polarization angle within if S [0,27r], and the (irrelevant) 
initial phase within ^g S [0, 27r]. 

The sky-position for all signals is fixed (without loss 
of generality) to that of Scorpius X-1, namely {a, S) = 
(4.276, —0.273) rad, and we assume gapless data from the 
LIGO Hanford detector (H) [3]. 

The random signal phase-evolution parameters Xg € V 
are generated by drawing them from uniform distribu¬ 
tions over the ranges: 


/ e [50, 1000] Hz , 
a sinj 


Q-n — 


G [1, 5]s, 


PG[Po- dP, Po + dP ], 


tp G 


P P 

^mid 5 ^mid T 


(79) 


£(mo,m) 


Mo - M 
0.5 (fj-o + m) ' 


(78) 


logige e [-5, lOgig(0.9)] , 
w e [0, 27r] rad. 


which is bounded within e G [—2,2]. For small values 
jej <C 1 this agrees with the asymmetric definitions of 
relative errors, such as eg = (/rg — m)/Mo, with domain 
eg G (—oo, 1], and which is related to the symmetric rel¬ 
ative difference e as eg = . 


A. Monte-Carlo software-injection method 

The general algorithm used for these software-injection 
tests is the following: 

1. Pick random signal amplitude parameters A and 
phase-evolution parameters As from suitable priors. 

2. Generate a phase-parameter offset (5A = At — Ag by 
finding the closest lattice template At to the signal 
point As. 

3. Gompute the phase metric gij{Xs) and metric mis¬ 
match /i(As, At) = gij SX'‘SX^ according to Eq. Q. 

4. Generate a (noise-free) data-set containing the sig¬ 
nal X = h{t;A,Xs)- Gompute the (coherent or 
semi-coherent) J^-statistic at the injection point, 
P{x;Xs), and at the template, J^(a:;At). In noise- 
free data we have J'(x;A) = E[A'(x;A)], therefore 
we can obtain p^(.A, As;A) via Eq. 0 01' Eq- 0, 
respectively, which yields the measured mismatch 
/ig via Eq. 0. Injection and recovery is performed 
using the LALSuite software package [49] . 

This procedure is applied to test each of the four types 
of metric detailed in Secs. ITV^ and lTVC] 


The orbital period “scale” Pg is fixed for each set of soft¬ 
ware injections (specified later), while the corresponding 
sampling range dP = Pg dH/(27r) is given in terms of 
a range dCl in orbital angular velocity, which is chosen 
as dn ~ 1/-y/gno, corresponding to a roughly unity mis¬ 
match along H. The motivation for this construction is 
twofold: on one hand we need to control the scale of the 
orbital period in order to ensure the appropriate short- 
segment (AT ^ P) or long-segment (AT ^ P) limit 
is satisfied in the tests. On the other hand, we wish to 
randomize the period over a range larger than the typi¬ 
cal template-bank spacings in order to fully sample the 
Wigner-Seitz cell of the template bank. 

For the metric tests presented in this section we employ 
a Demodulation method of computing the P-statistic, us¬ 
ing Short Fourier Transforms (SETs) of length Tgft = 5 s. 
Due to the linear phase model over each SET employed 
by the search code, this corresponds to a maximum er¬ 
ror in phase of JAc/ij ^ 0.07 rad over the investigated 
parameter space (see Appendix [Cj for more details). 


Generating template-bank offsets SX 

Sampling suitable offsets i5A for metric testing has 
proved to be a subtle and difficult point in previous met¬ 
ric studies [ini [3D] . Here we employ a recent innovation 
for generating “natural” offsets from a virtual template 
bank (e.g. see [3T| for more details). This method can be 
applied whenever the metric gij is constant (i.e. indepen¬ 
dent of A), which allows for constructing a lattice-based 
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template bank for a given maximal mismatch /Xmax of 
Eq. (12), namely 

Mmax > 5A* <7,, Sik , (80) 

where 5j3^ = i5A* Ai in terms of the Cholesky decomposi¬ 
tion of the metric, i.e. gij = Ai Sik Aj^. Further rescaling 
(5g* = (TZ/y/fimax) 5(3^ yields the corresponding Euclidean 
covering equation 

Sq" Sij Sq^ < , (81) 

in terms of the covering radius TZ. This requirement can 
be satisfied by lattices [311 EO] such as, for example, the 
simple n-dimensional hyper-cubic lattice Z„ (with cover¬ 
ing radius of TZz„ = ^/nl2), or the highly efficient cov¬ 
ering lattice A'^. For both Z„ and A* lattices, efficient 
algorithms exist for finding the closest lattice point qt, 
satisfying Eq. (811 for any given point qs € K", with 


Sq = qt — qs- For example, for Z„ this is trivially given 
by component-wise rounding, i.e. qt = round[gs]- We can 
use this construction to: (i) easily find the closest lattice 
template At to any given As, (ii) transform into lattice 
coordinates, i.e. qs = q{Xs), (hi) find the closest lattice 
point qt and invert back, i.e. At = X{qt)- The resulting 
distribution of offsets (5A = At — As uniformly samples 
the Wigner-Seitz cell of the corresponding lattice, which 
is the natural offset distribution for uniform signal prob¬ 
ability over the parameter space Xs G V. 

Note that while the coherent short-segment metric 
of Eq. (57) is explicitly flat (i.e. constant), this is not true 
for the other three cases, namely the semi-coherent short- 
segment metric of Eq. (62), the coherent long-segment 


metric of Eq. (47), and the semi-coherent long-segment 


metric of (48). However, it is easy to see that diagonal¬ 


rescaling via Eq. (40) achieves explicit flatness in all these 
three cases, namely by working in terms of rescaled co¬ 
ordinates 


SX^ = ^tSX\ 


(82) 


It might seem that the terms involving Ama = tmid ~ 
tasc depend on the coordinate taso but this effect can be 
neglected in all these cases, as we always assume Tobs 
P, and by choosing a gauge on tp as in Eq. (79), we can 


effectively approximate Ama/Tobs ~ 0 in the metric. 

The following metric tests use a Zg lattice and a max¬ 
imum mismatch of /tmax = /^max = 0.3. This lattice is 
chosen purely for simplicity, and only serves to generate 
realistic signal offsets for testing the metric. Using a dif¬ 
ferent lattice (e.g. ) would slightly change the distribu¬ 

tion of sampled offsets, but would be inconsequential for 
the metric-testing results. The maximal-mismatch value 
of Aimax = 0.3 represents a fairly “typical” value for real¬ 
istic searches, but is still small enough so that measured 
P-statistic losses are only minimally affected by higher- 
order corrections compared to the metric approximation 
(e.g. see |29l |30]). Here we are mostly interested in the 
accuracy of the metric within its range of applicability, 
while nonlinear deviations from the metric approxima¬ 
tion would warrant a separate study. 


B. Results in the long-segment regime 


1. Coherent long-segment metric g^^ 


For these tests we fix the “scale” of the orbital pe¬ 
riod to Pq = 19 h (similar to Scorpius X-1), and vary 
the segment length AT from 3 days to 39 days in steps 
of 4 days. We choose the parameter space range for 
the orbital velocity H as two times the largest possible 
metric spacing occurring over the parameter space, i.e., 
dCl = 2 max[(g^Q)“^/^] ~ 1.2 x 10“^ s“^. The resulting 
parameter-space half-width in Period is dP ~ 90s, and 
therefore we have at least AT/P >3.8. 

The results of these tests are shown in Fig. for the 
total 20 000 trials performed (2 000 trials for each AT 
value). We note that the performance of the parameter- 
space metric starts to degrade below AT/P < 10 for 
low-eccentricity orbits [panel (6)], and becomes generally 
poor for high-eccentricity orbits above e > 10“^ [pan- 
ela (a) and (c)]. Panels (d) and (e) show the agreement 
between the measured mismatch distribution and the ex¬ 
pected Zg-lattice distribution. In panel (/) we see that 
the phase-metric approximation of Eq. ^ agrees bet¬ 
ter for small mismatch values and develops a slight ten¬ 
dency to over-estimates the actual loss Jlq for higher mis¬ 
matches. This is a general feature seen in all four cases 
tested here, and qualitatively agrees with a similar effect 
seen in tests of the all-sky metric for isolated CW signals 

Hi 130]. 


2. Semi-coherent long-segment metric 


The scale for the orbital period used here is Pq = 2 h, 
with segments of fixed length of AT — 1 day and varying 
Tobs- In total we used 10 values of T^-ps'- Iday, 10 days, 
and 30 days up to 100 days in steps of 10 days. 

The parameter-space extent in orbital velocity is the 
same as considered in the previous section, i.e., dXl ^ 
1.2 X 10“^ s“^, resulting in a parameters-space half-width 
for the period of dP ^ Is. Therefore AT/P > 12 is 
satisfied for all trials. 

The results of these tests are shown in Fig. for the 
total of 20 000 trials performed (2 000 trials for each Tobs 
value). From panels (a, 6, c) we see that the agreement 
between measurements and predictions is generally good 
in the low-eccentricity regime e < 10“^, while rapidly 
degrading at higher eccentricites. The case Tobs = AT = 
1 day corresponds to the single-segment coherent case, 
and antenna-pattern effects are expected to play a role 
for observation times of order a day. 
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FIG. 1. Coherent long-segment regime: results for metric tests on of Sec. 
varying segment length AT from 3 to 39 days in steps of 4 days, (a) Relative error £{Jlo, m) versus eccentricity e. The dashed 
vertical line denotes the boundary between a “low-e” range e G [10“®, 0.1] (left-column plots) and a “high-e” range e G [0.1, 0.9] 
(right-column plots). (b,c): Relative error e vs observation time AT/P. The solid lines in panels (a, b, c) denote the median 
value, the error bars correspond to the 25th-75th percentiles, the circles and the stars denote the 2.5th and 97.5th percentiles, 
respectively. (d,e): Mismatch histogram of measured mismatches /io//lmax and theoretical distribution in a Ze lattice. (f,g): 
Measured mismatch Jlq versus predicted Jl. 
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FIG. 2. Semi-coherent long-segment regime: results for metric tests on ggj. 



for orbital period scale of Po = 2 h 


and fixed segment length of AT = Iday. (a) Relative error e(/lo, /I) versus eccentricity e. The dashed vertical line denotes the 
boundary between a “low-e” range e G [IQ-®, 0.1] (left-column plots) and a “high-e” range e G [0.1, 0.9] (right-column plots). 
(b,c): Relative error e vs observation time Fobs. The solid lines in panels (a,b,c) denote the median value, the error bars 
correspond to the 25th-75th percentiles, the circles and the stars denote the 2.5th and 97.5th percentiles, respectively. (d,e): 
Mismatch histogram of measured mismatches /io/)lmax and theoretical distribution in a Zg lattice. (f,g): Measured mismatch 
'po versus predicted p. 
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C. Results in the short-segment regime 

1. Coherent short-segment metric 


There is an important technical difference in the injec¬ 
tion algorithm in this case, as the metric is expressed in 
M-coordinates of Eq. (58) instead of physical coordinates 
A. Here the physical coordinates of an injected signal As 
are converted into it-coordinates, rts,fc = Uk{Xs), and then 
the closest lattice template Ut^k is found. In principle one 
could try to convert this back into physical coordinates 
using the expressions given in Appendix]^ but this is not 
possible in all cases, as some templates in u-coordinates 
do not correspond to physical coordinates and therefore 
cannot be inverted. In order to circumvent this problem, 
we use the fact that the w-coordinates essentially describe 
an isolated CW signal with n — 1 = 5 spindown values. 
Hence, the perfectly-matched J^-statistic is computed as 
usual at J^s = •A(As), but the J^-statistic for the mis¬ 
matched template is computed at the isolated spindown 
location defined by the u-coordinates, namely by setting 


/f) = ut,k+i for fc = 0,... 5, i.e. Tt = 




The orbital period scale is set as Pq = 80 days, and 
the resolution for the orbital angular frequency is taken 
to be dH = 2 X 10“^ s“^, and hence P S [62, 97] days. 
We sample segment lengths AT ranging from 2 days up 
to 25 days in steps of half a day, exploring the range in 
relative segment length of AT/P ^ [0.02, 0.4], using a 
total of 94 000 trials performed (2 000 trials for each AT 
value). 

In order to better understand these results, it is in¬ 
teresting to consider the effective template dimension for 
the different software injections, which we can quantify 
in terms of the highest nonzero component Ut.fc' of the 
closest template found for an injection, namely 


k' = max k with [Mt.fcj > 0 . (83) 

This quantity is plotted Fig. for the injections per- 



AT/P 


FIG. 3. Effective template dimension k' of Eq. (831 versus 
segment length AT /P for the software-injection studies shown 
in Fig. Solid line: mean value, error bars: maximum and 
minimum of k'. 


formed for Fig. giving the range of effective template 
dimensions k' as a function of relative segment length 
AT/P. As expected, we see that the effective template 
dimension increases with AT/P and eventually would 
require more than 6 M-coordinates for AT/P > 0.2, cor¬ 
responding to the region where the u-metric performance 
is seen in Fig. |^b) to start to degrade. We also see that 
the average template-bank dimension is not constant and 
generally less than 6, which explains the discrepancy be¬ 
tween the sampled mismatch distribution in Fig. |^d) 
compared to the expected Zg-lattice distribution. 

The results of the injection tests are shown in Fig. 
From panels (6) and (c) we see that the approximation in 
terms of 6 rt-coordinates used here starts to break down 
in the range AT/P > [0.25,0.35], as anticipated from 
Eq. (60) and Fig. Also, contrary to the other three 
limits tested, the small-eccentricity approximation only 
seems to work well for eccentricities up to about e 3 x 
10“^, after which it rapidly deteriorates and increasingly 
loses predictive power beyond e > 10“^ [see panel (a)]. 


2. Semi-coherent short-segment metric 

The scale for the orbital period here is fixed to Pq = 
10 days, and we used a fixed segment length of AT = 

1 day. The parameter space extent in terms of Q is dfl = 

2 max[(gQQ)~^^^] 6-6 x 10“® s“^ resulting in dP ~ 

2.2 hours and therefore we have at most AT/P < 0.12. 
We use a varying Tobs ranging from 30 up to 100 days in 
steps of 10 days. The results of these tests are shown in 
Fig. for the total 16 000 trials performed (2 000 trials 
for each Tobs value). We note that the validity extent of 
the semi-coherent metric is constantly acceptable for low- 
eccentricity orbits [panel (6)] and degrades (albeit much 
less than in the other three cases) for higher eccentricity 
e > 10“^ [panels (a) and (c)]. Panels (d) and (e) compare 
the mismatch histogram of measured mismatches with 
the theoretical distribution in a Zg lattice. Furthermore, 
we also note that in this regime the metric mismatch 
/t tends to generally over-estimate /to somewhat more 
strongly than in the other cases [see panels (/) and (g)]. 


D. General discussion of metric testing results 


Summarizing the results presented in Figs. m we ob¬ 
serve that the agreement between metric mismatch pre¬ 
dictions and the measured relative loss in T-statistic is 
generally very good (typically no worse than ^ 10—30%) 
within its range of applicability. In the following we 
further quantify where various approximations start to 
break down. 

The first-order small-eccentricity approximation (see 
Sec. HIB I starts to noticeably degrade only above e > 0.1 
[Fig. [3a)> Fig. [^a) and Fig. [5];a)] except in the coher¬ 
ent short-segment r^ime where it breaks down above 
e ^ 5 X 10“^ [Fig. 0a)]. Further investigation of this 
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FIG. 4. Coherent short-segment regime: results for metric tests on of Sec. 


IV Cl 


using M-coordinates, for orbital period 


scale of Po = 80 days and coherent observation times in the range AT € [2, 25] days, (a) Relative error s(jlo, Ji) versus 
eccentricity e. The dashed vertical line denotes the boundary between a “low-e” range e G [10“®, 10“^] (left-column plots) and 
a “high-e” range e G [10“®, 10~^] (right-column plots). (b,c): Relative error e vs observation time AT/P. The solid lines in 
panels (a, b, c) denote the median value, the error bars correspond to the 25th-75th percentiles, the circles and the stars denote 
the 2.5th and 97.5th percentiles, respectively. (d,e): Mismatch histogram of measured mismatches )Io//imax and theoretical 
distribution in a Ze lattice. (f,g): Measured mismatch JIq versus predicted Jl. Note that in all plots except (b,c) the range of 
relative segment length was restricted to AT/P < 0.2, where the short-segment limit approximation is seen to be valid. 


















































































16 



eccentricity e 









5 ®® 


of Sec. 


IVC2 


for orbital period scale of 
(a) Relative error eifio, fl) versus eccentricity e. The dashed 
[10“®,0.1] (left-column plots) and a “high-e” range e € [0.1, 0.9] 


FIG. 5. Semi-coherent short-segment regime: results for metric tests 
Po = 10 days, and a fixed segment length of AT = 1 day. 
vertical line denotes the boundary between a “low-e” range e G 
(right-column plots). (b,c): Relative error e vs observation time Tabs- The solid lines in panels (a, b, c) denote the median 
value, the error bars correspond to the 25th-75th percentiles, the circles and the stars denote the 2.5th and 97.5th percentiles, 
respectively. (d,e): Mismatch histogram of measured mismatches /lo//lmax and theoretical distribution in a Ze lattice. (f,g): 
Measured mismatch /lo versus predicted 
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behaviour might be interesting, but is beyond the scope 
of this study. 

The (phase-) metric approximation Eq. ®> is expected 
from previous studies [29l [30] to be quite accurate (in 
suitable coordinates) up to values of mismatch of /r < 0.2, 
after which higher-order terms start to become more no¬ 
ticeable, which tend to reduce the actual measured mis¬ 
matches compared to the metric predictions. While the 
nonlinear regime is hardly explored here, this general 
trend can still be seen to some extent in panels (f) of 

Figs. [Dll Hi 

VII. SCORPIUS X-1 SENSITIVITY ESTIMATE 

We can use the metric expressions and template counts 
derived in Sec. to estimate the optimal achievable sen¬ 
sitivity of a semi-coherent search directed at Scorpius X- 

1 . 



FIG. 6. Sensitivity depth of Eq. ( |84| ) for Scorpius X-1 

torque-balance level Eq. (851 versus CW frequency /. Dif¬ 
ferent line styles correspond to different progressive aLlGO 
sensitivities (“Early”, “Mid”, “Late” and “Final”) |53H55| . 
For each detector sensitivity three different emission models 
are shown, namely top: “mountain” (/ = 2u), mid: r-mode 
(/ ~ 4i//3), bottom: precession (/ ~ u), respectively. 


A. Torque-balance level 


In order to quantify the sensitivity of a search inde¬ 
pendent of the detector noise floor, we can define the 
sensitivity depth ED as 


B. Best theoretically achievable sensitivity 


Short-segment regime 


jjC — 

Pfa ~ hC ’ 


(84) 


in terms of the harmonic mean over detector noise floors 
Sn, and the strain sensitivity (or upper limit) of 
a search at a certain confidence level C (i.e. detection 
probability) and false-alarm threshold pfa. This quantity 
(with dimensions of is useful as a simple and in¬ 

tuitive measure for how far below the noise floor a given 
search setup can reach. 

An interesting astrophysical model postulating torque 
balance between the accretion and CW emission mmm 
yields a predicted CW amplitude (assuming a NS with 
10km radius and a mass of 1.4M©) for Scorpius X-1 of 


Let us first consider the short-segment regime, as¬ 
suming a “best-case” search with segments of length 
AT = 4.7h (i.e. roughly 1/4 of the Scorpius X-1 orbital 
period, quoted in Table |l]), total timespan Tobs = 365d 
(therefore N = 1864 segments), 2 detectors [LIGO H 
and Livingston (L)], and a very fine search grid of av¬ 
erage mismatch /i = 0.01. Using the method outlined 
in |56] and implemented in EZ], we obtain a resulting 
sensitivity depth of ^ 56Hz“^/^. Assuming 3 

equal-sensitivity detectors [LIGO H, L and Virgo (V)], 
this value increases to ^ 68Hz“^/^. 


Long-segment regime 


3.5 X 10 


-26 


300 Hz 


(85) 


where iz is the (unknown) NS spin frequency. 

Figure]^ shows the minimum required sensitivity depth 
pScoXi ].gach the Scorpius X-1 torque-balance limit 
of Eq. ( |8^ for different emission models (“mountain” 
deformation with / = 2 z/, r-mode emission with / « 
4z//3, and precession with / « v), assuming different 
aLIGO sensitivity curves (“Early”, “Mid”, “Late” and 
“Final”) [53H55]. 

For all the following sensitivity estimates, we will as¬ 
sume a detector duty cycle of 80% (following [55]), a 90% 
confidence level (i.e. detection probability), and a false- 
alarm level of Pfa = 10“^° (appropriate for a first-stage 
wide-parameter search). 


In the long-segment regime of Sec. IV B[ the maxi¬ 
mal segment length would be restricted by the astro- 
physical concern of “spin wandering”, namely a stochas¬ 
tic variability of the spin frequency due to variations 
in the accretion rate. There is substantial astrophys¬ 
ical uncertainty [13 dHl EEj about the details of this 
process and its magnitude, which is beyond the scope 
of this work. We consider assumptions roughly sim¬ 
ilar to those given in the recent Scorpius X-1 MDC 
which posited a random frequency derivative of or¬ 
der I/I < fs ^ 10“^^ Hz/s, changing on a timescale of 
tg ^ 10® s. Over this timescale, the frequency drift would 
therefore be 6fs < fgts 10“® Hz. Over a total ob¬ 
servation time Tobs, this can be modelled as a random 
walk with expected total drift (with respect to the mid- 
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TABLE I. Assumed Scorpius X-1 system parameters and un¬ 
certainties. Ranges of the form A ± AA denote the mean and 
Gaussian 1-sigma uncertainty, while [Ao,Ai] denotes a uni¬ 
form probability range. The time of ascending node tasc was 
computed from the originally-measured time of inferior con¬ 
junction To [48] via tasc = To —P/4. The projected semi-major 
axis Op is related to the projected radial maximal velocity Ki, 
appearing in |^, by flp = Ki P/{2 -kc). 


Parameter 

Value 

Ref. 

Op (s) 

1.44 ±0.18 

[Hisn] 

tasc (GPS S) 

897753994 ± 100 

mi 

P(s) 

68023.70496 ± 0.0432 

m 

e 

0) or (0.033 ±0.018) 

m 

uj (rad) 

[0, 27r] 

m 


point in time) of order Sfs •\/"^obs/(2ts) Sfs- The maxi¬ 
mal coherent segment length could therefore be estimated 
roughly by the requirement that the total frequency drift 
should be less than the frequency resolution ^ 1/AT of 
the search, in order to avoid any significant loss of SNR. 
This yields the constraint 

AT<V2 (^Tobs ts fs) . (86) 

Assuming Tobs = lyear and the above MDC spin¬ 
wandering model, we find AT < 3 d. Given the sub¬ 
stantial uncertainty in these parameters, for comparison 
we also consider a more optimistic scenario of AT < 10 d. 

The best achievable sensitivity for AT = 3d and 
N = 120 segments, assuming a small mismatch of 
/r = 0.01, can be estimated for a 2-detector network 
as ~ 105and for 3 detectors we find 

T^le-io ~ 127 Hz" 

The best sensitivity that can be achieved for AT = 
10 d and N = 36 segments, and fj, = 0.01, can be esti¬ 
mated for a 2-detector network as ^ 134Hz"^^^, 

and for 3 detectors we find T>®g?iQ ^ 163Hz"^/^. 

In the following we will focus only on the long-segment 
regime, and assess the required computing power as well 
as the resulting sensitivity depth that can actually be 
reached in practice. 


C. Scorpius X-1 search parameter space 


We assume the Scorpius X-1 parameters and uncer¬ 
tainties given in Table |T] Note that, contrary to the Scor¬ 
pius X-1 MDC Hn] and previous searches [111, ^Iso 
allow for a nonzero uncertainty on the eccentricity, in 
addition to the circular-orbit assumption that is also not 
ruled out by current observations [59j . 

Assuming 3-sigma ranges for the Scorpius X-1 pa¬ 
rameter uncertainties, the resulting per-dimension tem¬ 
plate numbers A//,, at maximal mismatch Umax = 0-1 of 
Eq. (70) are shown in Fig. We see that, for the poten- 



0 200 400 600 800 1000 1200 1400 

CW Frequency / (Hz) 

FIG. 7. Effective number of (coher ent or semi-coherent) tem¬ 
plates per dimension A// of Eq. ( |70| ) (in the long-segment 
regime) as function of GW frequency / for the Scorpius X-1 
parameter space given in Table]!] using ±3 sigma ranges. The 
template numbers are evaluated for ^max = 0.1 and mean val¬ 
ues of ftp and fl. Note that assuming a 27r uncertainty in cu 
entails = A/)?7r. The template count jVn is the only 

one depending on the observation time T and is shown for 
T= Id, lOd, 90d. 


tial uncertainty on eccentricity e of Table [Tj one would 
have to include e and w as effective search dimensions 
in the long-segment limit. On the other hand, the as¬ 
sumption of full circularization e « 0 is not ruled out by 
the current observations and therefore also constitutes a 
reasonable alternative model to investigate. 

The only dimension where the resolution depends on 
the timespan T is H. We see, for example, that for co¬ 
herent times AT < 10 d the uncertainty in H will typi¬ 
cally not be resolved by the coherent metric (as long as 
Tmax > 0.1), as A/n(T < lOd) < 1. In this case one 
could simply exclude this search dimension in the coher¬ 
ent template count A/ of Eq. (13). On the other hand, 
the semi-coherent fine grid (assuming Tobs > 90 d and 
Tmax < 0.1) can typically resolve this dimension over 
some frequency range where Afn{T > 90d) > 1. In both 
cases, however, this depends on the resulting mismatch 
parameter fj-max, which is itself an optimization param¬ 
eter. Therefore, as seen in Fig. one potentially has 
to deal with situations where the D-dimension will be 
unresolved at lower frequencies but resolved at higher 
frequencies. One way to obtain a correct template esti¬ 
mate in such a case would be to determine if a transition 
frequency exists, where the H-dimension changes from 
unresolved to resolved, then stitch together the respec¬ 
tive correct template counts over each frequency range. 
Alternatively, one can simply split the frequency range 
into narrow slices, such that the dimensionality can be 
assumed constant over each slice, and sum up the respec¬ 
tive number of templates. 

Note that due to the large separation Ama between the 
orbital epoch Gsc in Table [l] and the aLIGO gravitational 
wave search epoch, we have to consider the potential in¬ 
crease in uncertainty on tasc, as discussed in [IS]. As dis- 
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cussed in Sec. V A this is not a concern for a grid where 
both and tasc are fully resolved by the metric, as is 
expected to be the case of the semi-coherent grid. On 
the other hand, the coherent grids will typically not re¬ 
solve the O dimension, and therefore an increase of about 
2.5-3 in uncertainty on tasc would have to be taken into 
account. However, for the present study, we assume (as 
seems likely [IH]) that there will be further observations 
on Scorpius X-1 closer to the aLIGO epoch, which will 
again constrain the uncertainty on tasc to a level similar 
(or better) than that given in Table |l] 


D. Computing cost model 


As described in 1^. a semi-coherent StackSlide search 
(introduced in Sec. IIBI for a fixed parameter space has 
the following tuneable parameters: the segment length 
AT, the number of segments N, the maximal template- 
bank mismatch for the (semi-coherent) fine grid )Imax, 
and for the per-segment (coherent) coarse grids ptmax- 
The aim is therefore to maximize the resulting sensitivity 
over this space of search parameters, under the constraint 
of fixed total computing cost Ctot(ltmax, Mmax, A^, AT) = 
Co. 

The total computing cost Ctot of the StackSlide T- 
statistic can be written as EDI 


Ctot(Mmax, /imax, A^, AT) = C -|- C , (87) 


where C(/i„iax, A^, AT) is the cost of N per-segment coher¬ 
ent T-statistic searches over a coarse grid with maximal 
mismatch Jlraax, and C{'flyna,x, N, AT) is the cost of in¬ 
coherently summing (and interpolating) these T)-values 
across all segments on a fine grid with maximal mismatch 

/^max’ ^ 

The coherent cost for the T-statistic computation can 
be expressed as 

C = A A/'OIn.ax, at) Adet ?! (AT), (88) 


where Adet is the number of detectors, and ci is the (per- 
detector) computing cost per template, which depends on 
the algorithm used to compute the T-statistic. For the 
(generally slower) SFT-based demodulation method [61] . 
this can be expressed (for gapless data) as 


ef\AT)=^, 


AT 


(89) 


where « 4 x 10“® s is an implementation- and 
hardware-dependent fundamental computing cost for 
binary-CW templates (measured by timing the current 
T-statistic LALSuite code |49]1. 

Using instead the more efficient Fast Fourier Transform 
(FFT) based “resampling” method (DHUSl], we find the 
time per template ci to be approximately constant (for 
search frequency bands of > 10^ frequency bins) as 

« 3 X IQ-^s. (90) 


The semi-coherent cost of summing across segments can 
be expressed as 

C = AA(Ai„,ax,AT,A)do, (91) 

where Cq « 5 x 10“® s is an implementation- and 
hardware-dependent fundamental computing cost of 
adding one value of 2 Tt, for one fine-grid point (measured 
by timing the LALSuite HierarchSearchGCT code [49] 1. 

Note that the sensitivity estimates presented in Ta¬ 
bles |l^|V^ are using “Einstein@Home months” (EM) as 
a computing-cost unit. This unit is defined as 12 000 core- 
month^ on a CPU that achieves the above-quoted tim¬ 
ings Co and ci. Empirically this corresponds to the per¬ 
formance of the “average CPU” currently participating 
in Einstein@Home, and is comparable to roughly an Intel 
Core i7-2620M or Intel Xeon E3-1220 v3. 


E. Sensitivity estimates 


With the computing-cost model at hand we can em¬ 
ploy the semi-analytic machinery of [20] to find optimal 
StackSlide setups for given computing cost. We can then 
estimate the corresponding sensitivity of each setup with 
the fast and accurate method developed in [56] . Both the 
optimization and sensitivity-estimation algorithms have 
been implemented in m and are used here. 

Due to the rapidly increasing computing cost with fre¬ 
quency [see Eqs. (72)-(75)], it is important to limit the 
frequency range searched in some way. Hence, we fo¬ 
cus on constructing example setups that beat the torque- 
balance level (shown in Fig. over their respective fre¬ 
quency range. 

We consider the two different constraints on the seg¬ 
ment length AT due to spin wandering, as discussed 
in Sec. VH B[ namely AT < lOd and AT < 3d, re¬ 
spectively. For AT < lOd, the results for different T- 
computation methods and assumptions about eccentric¬ 
ity are summarized in TableO(for SFT-based demodula¬ 
tion and zero uncertainty on^ , Table HI (for FFT-based 
resampling and zero uncertainty on e), and Table IV 
(for FFT-based resampling and 3-sigma uncertainty of 
Ae = 0.087). 

Similarly, for AT < 3d, the results for different T- 
computation methods and assumptions about eccentric¬ 
ity are summarized in Table[V|(for SFT-based demodula¬ 
tion and zero uncertainty on^. Table VI (for FFT-based 
resampling and zero uncertainty on e), and Table VH 


(for FFT-based resampling and 3-sigma uncertainty of 
Ae = 0.087). 

In each table the first line shows a setup beating the 
f = 2v torque-balance depth for the “Final” aLICO sen¬ 
sitivity curve of Fig. [^ over the given frequency range. 


^ Counting only the 30%-50% CPUs currently devoted to gravita¬ 
tional wave searches, and including the double-computation of 
all results for validation. 
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while subsequent lines illustrate slight variations when 
changing only one of the constraints (computing cost 
Cq, observation time Tobs or detectors, respectively). In 
the last line of each table we also show a setup beat¬ 
ing the f = 2v torque-balance level depth over the given 
frequency range assuming a 6-month science run using 
the “Mid” aLIGO configuration, which is currently esti¬ 
mated [SS] to take place in about 2016-2017. 

Summarizing these results, we find that the SFT-based 
demodulation method can only reach the torque-balance 
level assuming we can neglect the uncertainty in eccen¬ 
tricity (Ae = 0), shown in Tables [TT| and [V| The FFT- 
based resampling method, in the case of Ae = 0, can beat 
the torque-balance level over a wider frequency range up 
to ~ 630 Hz (for AT < 10 d) or up to ^ 500 Hz (for 
AT < 3d), eve n at substantially reduced computing cost 
(cf. Tables HI and |VI[ ), and can still beat the torque- 
balance limit (albeit only up to a smaller frequency of 
~ 160 Hz — 200 Hz, depending on the constraint on AT) 
in the case of substantial uncertainty (Ae = 0.087) on 
eccentricity, as seen in Tables [IV| and [VH| 

At a lower-false alarm probability of pfa = 10“^'*, cor¬ 
responding to a level at which one might set upper limits 
without any further followup, the corresponding sensitiv¬ 
ity depth would be roughly 10% lower compared 

to the values of quoted in Tables 


VH 


Note that some of the optimal setups include maxi¬ 
mal mismatches Pmax of order unity, which would violate 
the metric approximation discussed in Sec. |H C| In order 
to fully quantify this effect, further Monte-Carlo tests 
on the higher-order metric deviations would be required, 
but qualitatively we know from previous results [e.g. see 
[iniisni, Figs. 0(f),i (f), and metric tests we performed 
at Pmax = 0.5 that are not shown here] that in this regime 
the metric approximation generally tends to overestimate 
the actual measured mismatches. Furthermore, the dis¬ 
tribution of sampled mismatches will be peaked around 
their mean value, which for lattices is (/i) « 0.5 Pmax, 
and will therefore not be affected as much by these devi¬ 
ations. 


F. Comparison to previous sensitivities 


In order to put these sensitivities in context, we 
consider some previously-achieved sensitivities for Scor- 
pius X-1 CW searches. The first coherent Scorpius X- 
1 search presented in m was computationally limited 
to using a total of ~ 9h of data collected by LIGO 
during the second science run with the Livingston and 
Hanford detectors. Such a search achieved upper lim¬ 
its corresponding to a sensitivity depth of approximately 
T>95% ~ 4Hz“^/^, which is roughly consistent with the 
theoretical expectation [56] of « 4.5 Hz“^/^ for a 

fully-coherent two-detector search. The recent search on 
the fifth LIGO science run, which uses a semi-coherent 
sideband method m, achieved a substantially improved 


sensitivity depth of T>95% 30Hz“^/^ over a frequency 

band of 50 — 550 Hz. Note that the software-injected sig¬ 
nals used in the recent Scorpius X-1 MDC nni are found 
(over a frequency range of 50 — 1500 Hz) in the depth 
range T> € [2.6, 104] Hz“^/^ for the “open” injections, 
and V € [2.4, 58] Hz“^/^ for the “closed” injections. 


VIII. CONCLUSIONS 

The work presented here consists of three main parts. 

In the first part we re-derived the (coherent and semi- 
coherent) metric expressions for the long-segment and 
short-segment regimes, assuming a low-eccentricity bi¬ 
nary orbit. We found that (as has also been noted ear¬ 
lier), the long-segment regime does require refinement 
in the orbital angular velocity H, contrary to an earlier 
result in the literature m- We have extended these find¬ 
ings to allow for general offsets between the orbital ref¬ 
erence time epoch tasc and the observation epoch, which 
explicitly show that there is generally no increase in com¬ 
puting cost with increasing time passed since the origi¬ 
nal orbital reference epoch if both H and Gsc are fully 
resolved by the metric. 

In the second part we subjected the analytic metric 
expressions to extensive Monte-Carlo tests comparing 
their predictions against measured SNR loss in software- 
injection studies, and found robust agreement within 
their range of applicability. 

In the final part of this work we used the metric 
template counts to semi-analytically estimate the opti¬ 
mal achievable StackSlide sensitivity of a directed search 
for Scorpius X-1. We found that the predicted torque- 
balance limit could be reachable for the first time with an 
Einstein@Home search using data from a 6-month “Mid” 
aLIGO run (currently planned for 2016-2017), and should 
be beatable over a substantial frequency range in “Fi¬ 
nal” aLIGO. However, the frequency range over which 
the torque-balance limit can be beaten depends strongly 
on the Scorpius X-1 parameter-space uncertainty, most 
notably the orbital eccentricity and the assumed charac¬ 
teristics of spin wandering. 

More work is required to better understand the mod¬ 
elling and effects of spin wandering, as well as the re¬ 
sulting constraints on search methods. It would seem 
desirable to develop a more robust statistic that prop¬ 
erly takes spin wandering into account, for example by 
marginalizing over the spindown uncertainty in each seg¬ 
ment before incoherently summing them. 
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TABLE II. Optimal search setups and corresponding sensitivity estimates for a Scorpius X-1 StackSlide search (with AT < 10 d) 
using the SFT-based demodulation algorithm, assuming negligible uncertainty Ae = 0 on the eccentricity. The SET length 
assumed is rsft(430 Hz) = 240 s and Taft(250 Hz) = 320 s (see Appendix [c| for /Xaft = 0.01). The hrst line shows a setup 
beating the “Final” / = 2iz torque-balance depth (see Fig.[^ over the given frequency range, and subsequent lines illustrate the 
results for slightly different constraints on computing cost Co, observation time Tobs and detectors (IFOs). The last line shows 
a setup beating the f = 2v torque-balance level depth over the given frequency range with a “Mid” aLIGO configuration. In 
all cases the template-bank dimensions are found asn — 3 and n = 4. The computing-cost unit “EM” corresponds to running 


on Einstein@Home for a month, and is 

explained in more 

detail in Sec. |VH D 





/ 

(Hz) 

Co 

(EM) 

IFOs 

Tobs 

(d) 

N 

AT 

(d) 

/^max 

/^max 

Af 

Af 

C/C 

-ri90% 

^le-lO 

(I/VHz) 

[20,430] 

12 

HL 

360.0 

43 

8.30 

0.71 

0.04 

3.4x10^^ 

7.5x10^® 

22 

100 

[20,430] 

6 

HL 

360.0 

83 

4.36 

0.74 

0.06 

1.7x10^^ 

2.4x10^® 

18 

86 

[20,430] 

12 

HL 

180.0 

18 

10.00 

0.50 

0.02 

6.9x10^® 

1.5x10^'^ 

27 

91 

[20,430] 

12 

HLV 

360.0 

62 

5.77 

0.73 

0.04 

2.3x10^® 

5.1x10^® 

22 

112 

6 months science run with “Mid” aLIGO: 

[40, 230] 

6 

HL 

180.0 

18 

10.00 

0.18 

0.010 

4.8x10^® 

7.8x10^® 

26 

102 


TABLE HI. Optimal search setups and corresponding sensitivity estimates for a Scorpius X-1 StackSlide search (with AT < 
10 d) using the FFT-based resampling algorithm, assuming negligible uncertainty Ae = 0 on the eccentricity. In all cases the 
template-bank dimensionalities are found as n = 3 and n = 4. The notation is the same as in Table [H] 


/ 

(Hz) 

Co 

(EM) 

IFOs 

Tobs 

(d) 

N 

AT 

(d) 

/^max 

/^max 

AT 

Af 

C/C 

^90% 

■^le-10 

(1/^) 

[20, 630] 

12 

HL 

360.0 

36 

10.00 

0.04 

0.03 

1.1x10^® 

00 

X 

H-' 

o 

1.6 

129 

[20, 630] 

6 

HL 

360.0 

36 

10.00 

0.06 

0.05 

5.5x10^® 

3.8x10®®" 

1.7 

127 

[20, 630] 

1 

HL 

360.0 

36 

10.00 

0.19 

0.12 

9.7x10^"^ 

5.6x10®® 

2.1 

117 

[20, 630] 

12 

HL 

180.0 

18 

10.00 

0.03 

0.02 

2.5x10®® 

1.4x10®® 

2.1 

105 

[20, 630] 

12 

HLV 

360.0 

36 

10.00 

0.05 

0.03 

7.5x10®® 

7.2x10®®" 

1.9 

155 


6 months science run with “Mid” aLIGO: 


[40, 250] 

6 

HL 

180.0 

18 

10.00 

0.006 

0.004 

1.1x10®® 

6.9x10®®" 

2 

107 

TABLE IV. Optimal search setups and corresponding sensitivity estimates for a Scorpius X-1 StackSlide search (with AT < 
10 d) using the FFT-based resampling method, assuming a pessimistic uncertainty Ae = 0.087 on eccentricity (see Table [I|. In 
all cases the template-bank dimensionalities are found as h = 5 and n — 6. The notation is the same as in Table [H] 

/ 

(Hz) 

Co 

(EM) 

IFOs 

Tobs 

(d) 

N 

AT 

(d) 

/^max 

/^max 

Af 

Af 

C/C 

^90% 

■^le-10 

(1/^) 

[20, 200] 

12 

HL 

360.0 

36 

10.00 

1.18 

0.54 

1.3x10®® 

5.7x10®"^ 

2.6 

60 

[20, 200] 

6 

HL 

360.0 

36 

10.00 

1.56 

0.68 

6.3x10®® 

2.8x10®"^ 

2.7 

55 

[20, 200] 

12 

HL 

180.0 

18 

10.00 

0.88 

0.35 

2.6x10®® 

10®® 

3 

56 

[20, 200] 

12 

HLV 

360.0 

36 

10.00 

1.37 

0.55 

8.6x10®® 

5.2x10®"^ 

3 

68 


6 months science run with “Mid” aLIGO: 


[60,140] 

6 

HL 

180.0 

18 

10.00 

0.57 

0.21 

1.3x10®® 

5x10®"^ 

3.1 

74 
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TABLE V. Optimal search setups and corresponding sensitivity estimates for a Scorpius X-1 StackSlide search (with AT < 3d) 
using the SFT-based demodulation algorithm, assuming negligible uncertainty Ae = 0 on the eccentricity. The SET length 
assumed is raft(420Hz) = 250s and rsft(150Hz) = 410s (see Appendix [c| for = 0.01). The first line shows a setup beating 
the “Final” f = 2u torque-balance depth (see Fig.[^ over the given frequency range, and subsequent lines illustrate the results 
for slightly different constraints on computing cost Co, observation time Tabs and detectors (IFOs). The last line shows a setup 
beating the f = 2v torque-balance level depth over the given frequency range with a “Mid” aLIGO configuration. In all cases 
the template-bank dimensions are found as n = 3 and n = 4. 


/ 

(Hz) 

Co 

(EM) 

IFOs 

Tobs 

(d) 

N 

AT 

(d) 

/^max 

/^max 

N 

N 

C/C 

^90% 

■^le-10 

(l/\/H^) 

[20,420] 

12 

HL 

360.0 

120 

3.00 

0.34 

0.03 

3.5x10^^ 

4.3x10^® 

14 

94 

[20,420] 

6 

HL 

360.0 

120 

3.00 

0.54 

0.05 

1.8x10^^ 

2x10^® 

15 

87 

[20,420] 

12 

HL 

180.0 

60 

3.00 

0.21 

0.02 

7.1x10^^ 

7.4x10^® 

16 

80 

[20,420] 

12 

HLV 

360.0 

120 

3.00 

0.44 

0.04 

2.4x10^® 

3.6x10^® 

16 

110 


6 months science run with “Mid” aLIGO: 


[50,150] 

6 

HL 

180.0 

60 

3.00 

0.03 

0.003 

5.7x10^® 

4.3x10^® 

14 

85 


TABLE VI. Optimal search setups and corresponding sensitivity estimates for a Scorpius X-1 StackSlide search (with AT < 3 d) 
using the FFT-based resampling algorithm, assuming negligible uncertainty Ae = 0 on the eccentricity. In all cases the 


template-bank dimensionalities are found 

as n = 

3 and n 

= 4. The notation 

is the same as 

in Table |II] 



/ 

Co 

IFOs 

Tobs 


N 

AT 

/^max 

/^max 

N 

Sf 

C/C 

'D90% 

■^le-10 

(Hz) 

(EM) 


(d) 



(d) 






(1/Vh^) 

[20, 500] 

12 

HL 

360.0 


120 

3.00 

0.02 

0.02 

3.2x10^® 

2.4x10^^ 

1.6 

102 

[20, 500] 

6 

HL 

360.0 


120 

3.00 

0.04 

0.03 

1.6x10^® 

l.lxlO^’’ 

1.7 

101 

[20, 500] 

1 

HL 

360.0 


120 

3.00 

0.12 

0.07 

2.9x10^"^ 

1.7x10^® 

2.1 

96 

[20, 500] 

12 

HL 

180.0 


60 

3.00 

0.01 

0.01 

6.7x10^® 

4.4x10^^ 

1.8 

84 

[20, 500] 

12 

HLV 

360.0 


120 

3.00 

0.03 

0.02 

2.3x10^® 

2.1x10^’’ 

1.9 

123 

6 months 

science run 

with “Mid 

’ aLIGO 









[50,190] 

6 

HL 

180.0 


60 

3.00 

0.003 

0.002 

3.4x10^® 

2.1x10^’’ 

2 

85 

TABLE VII. Optimal search setups and corresponding sensitivity estimates for a Scorpius X-1 StackSlide search (with AT < 
3d) using the FFT-based resampling method, assuming a pessimistic uncertainty Ae = 0.087 on eccentricity (see TableIn 

all cases the template-bank dimensionalities are 

found as 

n = 5 and 

n — 6. The notation is 

the same as 

in Table |H] 

/ 

Co 

IFOs 

Tobs 


N 

AT 

/^max 

^max 

JV 


C/C 

^90% 

■^le-10 

(Hz) 

(EM) 


(d) 



(d) 






(1/Vh^) 

[20,160] 

12 

HL 

360.0 


120 

3.00 

0.76 

0.34 

3.8x10^® 

1.7x10^'^ 

2.6 

58 

[20,160] 

6 

HL 

360.0 


120 

3.00 

1.00 

0.44 

1.9x10^® 

8.3x10^® 

2.7 

50 

[20,160] 

12 

HL 

180.0 


60 

3.00 

0.57 

0.22 

7.8x10^® 

3.1x10^'^ 

3 

59 

[20,160] 

12 

HLV 

360.0 


120 

3.00 

0.88 

0.35 

2.6x10^® 

1.6x10^'^ 

3 

65 

6 months 

science run 

with “Mid 

’ aLIGO 









[70,110] 

6 

HL 

180.0 


60 

3.00 

0.34 

0.13 

3.9x10^® 

1.5x10^'^ 

3.1 

71 


Appendix A: Inverting the u-coordinates back into 
physical coordinates 

Let us define re-scaled coordinates 


Wk = 


Uk 


(Al) 


and consider the set of Eqs. 


(581 for Wk 


Wk{X). 


1. General 6D elliptic case 

Let us take linear combinations of the equations for 
W 2 , ■ ■ - w^ to express 


sin 

cos V'm 
sin 2'0i„ 

cos r/’m 


1 

3/ap 


{W4 + 4 ^ 2 ) , 


{W 5 + Aws), 

3/ap 


— j ^ [ 2 k{w 2 + W 4 ) + r]{w3 

—^[k(u;3 + W5) - 2??(w2 
Qjp C 


(A2) 

(A3) 

W 5 )], 
(A4) 

wa)] , 
lA5i 
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where we used + rf. We now insert these into 

the equation for wq to obtain 

Wq + 5u>4 + 4^2 = 0 , (A6) 

which yields a quadratic equation for with solution 


= - 



(A7) 


which requires 25U4 — 16 u 2 M 6 > 0 for real-valued 
Note that once we know VL, then all Wk are known numer¬ 
ically as well. Inserting Eqs. (A2 1 -( A51 into the equation 
for wi brings us to find 

Using sin^ i/'m + cos^ = 1, we obtain 

3/ap = a/( iC5 -I- 4 w3)2 -I- {wi + 4 ^ 2 )^ , (A9) 

where we note that for constant-frequency signals, we 
numerically have f = f- Hence, given / from the previ¬ 
ous equation, we easily get Op. Proceeding similarly for 
sin^ 2^/^in -I- cos^ 2'!/)ni = 1, we find 


1 


e = 


12 /a, 


+ '^5)'^ + 4( w 2 + ^4)^ ■ 


(AlO) 


We can simply invert Eqs. (A2),(A3) to obtain and 
therefore tasc- Finally, solving the linear system of equa¬ 
tions W 2 ,W 3 for K,r], yields 


1 


77 = 


[2(w 2 -/apSin'0m)sin27/;n 

4/Op 

+ {W 3 - fttp cos 1pm) COS 2lpm] 

;j4— [{W3 - /Op cos 1pm) sin 2lpm 
4/Op 


2{fap sinipm - W 2 ) cos2ip„ 


(All) 


(A12) 


which, using Eq. (|A10) and Eq. ([33|) gives us to. 


which numerically determines all Wk ■ We can further see 
that 



fap = ^wl + wl, 

(A18) 

and 

tan iPm = —, 

W 3 

(A19) 

yielding tasc, 

and finally 



/ = fapflcosipm + Ul ■ 

(A20) 


Appendix B: Maximal Doppler shift due to orbital 
motion 


Sometimes it is important to estimate the maximal 
Doppler shift the intrinsic signal frequency of a binary 
CW signal can undergo due to orbital motion. From the 
phase model of Eq. (291 we see that the instantaneous 
Doppler shift is 


\/l — cos E cos uj — sin E sin uj 


d(p/dt ^ 


dR 

27rf 


cdt 


= OpH 


< OpH- 


1 — e cos E 

sinEsinwl -I- | cosEcoswl 
|1 — ecosEj 

where we used the fact that |a -I- 6 | < |a| 4 
•\/l — < 1. We further observe that 

I cosEcoswl -I- I sin if sin wl = 

max{ I cos(if -f w) |, | cos(if — w) |} < 1 , 

and 11 — e cos E\ > 1 — e to obtain 


dcp/dt 


27rf 


- 1 


< 


OpH 

1 — e 


. (Bl) 

| 6 | and 

(B2) 

(B3) 


2. Special 4D circular case 

In the special circular case, we can obtain a simpler 


solution by inverting the equations for {wi,... W 4 } and 
setting K = 77 = 0 , namely 

/ 

Wi = -- fttpCOSIpm, (A13) 

W 2 = fap sinipm, (A14) 

W 3 = fap cos 1 pm , (A15) 

W 4 = - fap sinipm ■ (AI 6 ) 


By summing equations for W 2 and W 4 , we find W 2 +W 4 = 
0 , which is a quadratic equation for B with solution 


Appendix C: Maximal SFT length 


By using the SFT-based demodulation method to com¬ 
pute the J^-statistic, the computing cost per template 
increases linearly with the number of SETs used [see 
Eq. ( 881 ]. On the other hand, the maximal length of the 
SFT is limited by the linear-phase approximation over 
the duration of each SFT. In order to minimize the com¬ 
puting cost of this method, we therefore want to choose 
the longest possible SFT duration [see Eq. (89)] with 
an acceptable error in the linear-phase approximation. 

In order to estimate the maximal phase-error of the 
linear-phase approximation over an SFTs, we can con¬ 
veniently re-use the short-segment limit expressions, and 
simply estimate the phase error as 


B = 



(A17) 




I I _ r, \'U'2\ I ,tsft \ '1 - 2 rT,2 

|w2| — 27r ^ ^ f ^ , 


Tkt 
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(Cl) 
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and the corresponding mismatch is given via Eq. (571 as 
Msft ~ 922’^'^2 — Turning this around to express 

the maximal Tsft for a given maximal mismatch /Xsft, we 
obtain 


7?ft(/) < , (C2) 

Trap / 

where the largest values of the parameter space being 
searched should be used for Op, n and /. This constraint 
on Tsft is illustrated in Fig. as a function of search 
frequency / for Scorpius X-1 parameters of Table We 



CW Frequency / (Hz) 


FIG. 8. Maximal allowed SFT length Taft of Eq. (C2| for a 
Scorpius X-1 demodulation 7^-statistic search as function of 
frequency /, for three different tolerated maximal mismatches 

/^sft■ 


see that for Scorpius X -1 this limit is more stringent than 
the analogous constraint coming purely from the Doppler 
effect due to the detector motion, which typically results 
in a choice of Tgft = 1800 s in searches for isolated CW 
sources. 


Appendix D: Optimal StackSlide solution for 
degenerate computing cost function 


In the case of the FFT-based resampling algorithm 
for the JF-statistic of Eq. (90), we encounter a degener¬ 
acy in the computing-cost function that had not been 
considered in the original optimization study of [20j . 
Namely, when the search dimension D is resolved in the 
semi-coherent template bank, but unresolved in the per- 
template coherent banks, then we see from Eqs. (72)-(75) 
that the number of templates (both for including eccen¬ 
tricity n = 5,n = 6 , and for the circular case n = 3,n = 
4), scale as oc AT and Af oc NAT, respectively, where 
we assumed the ideal gapless case with 7 = A. From 
the computing-cost expressions Eq. ([^ and Eq. (911 we 


see that therefore the computing-cost functions take the 
form C oc NAT oc Tobs and C oc (AAT)^ oc T^^bs- 
this case the power-law coefficients in the formalism of 
[3D] are rj = 6 = 1 and rj = 5 = 2 (the resulting N co¬ 
efficients at fixed Tobs are therefore e = e = 0 ), which 
corresponds to a degenerate case that has not been an¬ 
alyzed previously. In this case one can only find a so¬ 
lution by constraining Tobs- However, maximization of 
sensitivity at fixed computing cost would then result in 
A —>■ 1 , i.e. a fully coherent search (except if by increas¬ 
ing AT the template bank eventually starts to resolve 
D and therefore breaks the degeneracy). In cases such 
as the Scorpius X-1 search considered here, there is an 
astrophysically-motivated upper bound on the segment 
length AT < 10 d, and we therefore need to express the 
optimal solution for constraints on both Tobs and AT. 
With A and Tobs fixed, we can only optimize sensitivity 
over the respective template-bank mismatches /imax and 
/imax- This is achieved simply by solving Eq. (91) of [30] . 
i-e. (/fmax/n)A/imax/n) = C/C, together with C -bC = Cq 
for //max and /imax- 
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